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, Abstract 

, The 5D Cosmological General Relativity theory developed by Carmeli 

■ reproduces all of the results that have been successfully tested for Ein- 

^ | stein's 4D theory. However the Carmeli theory because of its fifth dimen- 

£f~^ . sion, the velocity of the expanding universe, predicts something different 

f~* > ; ' for the propagation of gravity waves on cosmological distance scales. This 

\£) t analysis indicates that gravitational radiation may not propagate as an 

unattenuated wave where effects of the Hubble expansion are felt. In 

""^ ' such cases the energy does not travel over such large length scales but is 

^ji evanescent and dissipated into the surrounding space as heat. 

5-H • Key words: cosmology, Carmeli, gravitational waves, 5 dimensions, expand- 
ing universe 

1 Introduction 

5-H 



In recent decades, the search for gravity waves has intensified with large high 
powered laser-based interferometic detectors coming on line. See LIGO ^3] 
and TAMA for example. These detectors have already reached sensitivities 
that should enable them to "see" well beyond the local galactic Group. On the 
other hand, the Hulse- Taylor binary |17j ring-down energy budget is a precise 
test of general relativity and a clear indication of the existence of gravitational 
radiation, and it seems that the first direct detection is just a matter of time. 

In standard General Relativity the expanding universe has no impact on the 
properties of gravitational waves, except the the well known effect of redshift. 
However, in Carmeli cosmology the expansion of the universe (or redshift of 
the gravitational wave) manifests as a fifth dimension and in this paper we 
calculate the effect and how this might impact on a possible direct detection. 
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2 Cosmological General Relativity 



In the late 1990s Moshe Carmeli proposed a new cosmology, Cosmological Gen- 
eral Relativity (CGR). 0211310113 It is a generally covariant theory and extends 
the number of dimensions of the universe by the addition of a new dimension 
- the radial velocity of the galaxies in the Hubble flow. The Hubble law is as- 
sumed as a fundamental axiom for the universe and the galaxies are distributed 
accordingly. 

As a result we have a 5D spacetimevelocity universe with two timelike and 
three spacelike coordinates in the metric. The signature is then — — — 
The universe is represented by a 5-dimcnsional Ricmannian manifold with a 
metric and a line element ds 2 — g^dx^dx v . This differs from general 
relativity in that here the x 4 — tv coordinate is more correctly velocitylike 
instead of timelike as is the case of x° = ct, where c is the speed of light, a 
universal constant and t is the time coordinate. In this theory x 4 = tv, where r 
is also a universal constant, the Hubble-Carmeli time constant. The other three 
coordinates x , k — 1, 2, 3, are spatial and spacelike, as in general relativity. 

It has been shown that all the results predicted by general relativity and ex- 
perimentally verified are also predicted by CGR 0] . However in [S] Carmeli dis- 
cussed the one consequence that was not exactly reproduced, and that was grav- 
ity waves in 5 dimensions. The new metric resulted in a redshift dependence with 
a more general wave equation incorporating 5 dimensions (ct, x , x 2 , x 3 , tv). 



2.1 Linearized general relativity 

As is the usual practice in a weak gravitational field we write the metric 

9fj.u = T)fj,v + h^v (1) 

where the metric rj^v is Minkowskian but extended here to 5D from the usual 
4D Minkowski metric but with signature (+ — — — +). Here rj^ is perturbed 
due to gravitating sources with <C 1. 

A useful tool is to define the trace-reversed as 

/V = h^ v - ^q^h (2) 

where h = rf^h a p is the trace of h^ v . Consequently 

h^u = h^v -^V»uh (3) 

where h = rj al3 h a p and h = —h. 

Then the linearized Einstein field equations to first order in yield 

O V = -2 K T M „ plus rfPh^p = (4) 

where Q is the D'Alembertian operator in 5D and may be expressed as 
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For conservation of energy and momentum, excluding gravity, it follows from 
that 

Tl a0 T^ p = 0. (6) 

From and it is clear that is a generalized wave equation that 
reduces to 

O V = (7) 

in vacuum. 

So the gravitational waves depend not only on space and time but also on 
the expansion velocity of the source in the Hubble flow. Here v is the fifth 
co-ordinate, which represents the velocity of the expansion of the space through 
which the wave passes. 

The solution of (jJJ is the sum of the solution of the homogeneous equation 
(|7|) . and a particular solution. The following is the special time independent 
retarded solution in the absence of source-less radiation. The contravariant 
form is 16, 



where the source mass is located at x' and the potential measured at x. To 
evaluate the integral in (JSJ) , provided the measurement point determined by the 
vector x is far away from the source, a Taylor expansion of l/|x — x'| about x' 
= is taken retaining only the first two terms, 

1 1 x k x lk 

i T\ ~ - + -IT" 9 

|x — x'| r r J 

where r 2 = x k x k . 

The integral in (JHJ is then written as 

- M „ = _GM_ _ Gjt nsn afi ( 
r 2 r A 

where the dipole term has been eliminated by choosing the origin of the coor- 
dinates to coincide with the center of the source mass. The first term in (|10|) 
involves the integral J T 00 (x')(i 3 a;' = M identified with the source mass and 
the second term J x /k T l0 (x.')d 3 x' — l/2e kln S n where S n is the spin angular 
momentum of the system. Here k,l,n = 1, 2, 3 for the spatial coordinates. 

2.2 Wave equation in curved spacevelocity 

Now considering the time dependence again we can write as 
1 B 2 \ - 1 fl 2 - 

where the substitution v/c — > z has been made. Here z is the redshift of the 
wave, and the substitution is valid where z =< 0.1, which is approximately 400 
Mpc. We assume it is approximately valid beyond that. Now the solution to 
is the sum of the solution to the homogeneous equation (|ll|l . which is the 
usual gravity wave solution in general relativity, and a particular solution of 
||TJ , which has a redshift dependent source term. 
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In fact, in CGR, because the Hubble law is assumed a priori, the expansion 
velocity v (or gravitational wave redshift z) is not independent of r and depends 
on the matter density of the universe. In fact, in the case of CGR, Ijlll) can be 
written as 

1 d 2 1 B 2 \ - 

^-V 2 + ^^U- = 0, (12) 



with 



when the chain rule is applied. 

Let us look for a plane wave solution of the form 




(13) 



hT = e^ v cosk a x a , (14) 

where the 3-space co-ordinates are (x , x 2 , x 3 ). Here x 1 and x 2 are orthogonal 
to the direction of propagation a; 3 from source to detector. Here e^" is a constant 
tensor and k a is a constant vector. Therefore we look for a wave propagating 
in the r direction which would have k a — (ui/c, 0, 0, k r ). 

This means we can retain only the r derivative in V 2 and effectively re- write 
(I12fl in spherical co-ordinates as 



Id 2 d 2 2d Id 2 
c 2 dt 2 dr 2 r dr c 2 t 2 dz 2 



W" v = 0. (15) 



Eqs l|12|) and l|15f) are only valid where the Hubble law applies. When it 
doesn't apply, that is, where dr/dz = in (|13fl . (|12|l becomes the normal wave 
equation for gravity waves in source free regions. 

However where the Hubble law is applicable, in flat (i.e. f2 = 1) spacevelocity 
dr/dz = ct, which is the Hubble law in the zero distance/zero gravity limit. In 
the general curved spacevelocity, the form of the derivative is given by [H] 



l + (l-fi)— (16) 



where f2 = p/p c is the mass/energy density at some epoch expressed as a fraction 
of the 'critical' density, p c = 3/8ttGt 2 . 

Substituting into JTSJ with (JT3J we get 

Id 2 2d f 2 d 2 »1\ Tllv n 



dt 2 r dr c 2 t 2 \ dr 2 dr 

which only has dependence on r and t. This is a new equation depends on the 
surrounding matter density Q. 

When the gravity wave is very distance from the source and when r 3> 
ct I V 1 — Q the second term of <|17[) is much smaller than the term in curly 
brackets, we assume the second term negligible. Therefore (|17|) can be approx- 
imated for large r as 

Id 2 l-n f 2 d 2 d) \ Tllv , , 

?s? + ^rp +p * r =a (18) 
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The solution of H18|) can be obtained by separation of variables assuming a 
solution of the form oc R(r)e l( - UJt+kz ' ) . Substituting the latter into <(T%}l yields 

u 2 n - 1 

"■T + -tor = 0- !9 



with k z w and i?(r) = air -1 for r ^> cr/y/1 — 0. More generally i?(r) = 
^a„r _n a polynomial expression with an index n > 0. Furthermore by taking 
a hint from the solution of the usual heterogeneous equation shown in lfTU|) R(r) 
is determined as 

R(r) = -— -O(i) 3 . (20) 
r r 

Equation (|19|) is a resonance condition. For this solution, which spans the 
whole extent of the Universe, the Universe acts like a resonant mode with a 
characteristic scale radius El of 



„2 



x /\Rl\ = \/\^\= „;' (21) 



and resonance frequency 



. = ^HE. (22) 
r 

For values of r = 4.2 x 10 17 s and = 0.02 the scale radius is Rq ~ 4.13 Gpc 
and the characteristic frequency is Llj/2tt « 3.66 x 10 -19 Hz. 

When r cr/\/l — £1 the second and fourth terms of l|17|l are much smaller 
than the third, and hence can be neglected. Therefore l|17|) can be approximated 
as 

\d 2 i - n 2 a 2 



C 2 ^2 ■ C 2 T 2 r 9r 2 J - - ( 23 ) 

Now for a spherically symmetric expanding universe Carmeli E] obtains 
the relation between redshift and distance of the emitting source, 



t_ _ sinh^yi - n) 

where q = ((1 + z) 2 — 1)/((1 + z) 2 + 1). Assuming this relation also holds for 
gravity waves 1|18|) becomes 

I + sinh 2 (,VT^)|^) = 0, (25) 

This is now a wave equation with a solution of the form hP" oc e ll , k rr+ut) ^ 
which results in the following dispersion relation 

2 a; 2 /c 2 f26) 

* r sinh 2 (,yi— W ^ 

which can be approximated for z <C 1 as 
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Figure 1: The decay constant ck/lu as function of redshift, z shown from Eq. 
(|26|l as the solid curve 1 and from the approximated Eq. I|27[) as the broken 
curve 2. |ck/w| asymptotes to unity in the limit as v — > c or when z — > oo. 
The values of fi changes from fi < 1 to fi > 1 at some value of z but the sinh 
function is continuous changing into a sine function as the argument of the sinh 
function changes from real to imaginary 



When fi > 1 the wave number is real and approximately k r = lu/(cz\/Q — 1) 
and hence gravity waves propagate yet are dependent on redshift. When fi < 1 
the wave number is imaginary and the amplitude is attenuated with a decay 
constant k = ik r = uj /{cz\/l — fi). In Fig. ^ we plot the decay constant n 
normalized by uj/c, using k r from (123). From the fi gure it is apparent that the 
approximation of (|27|l is good for z < 0.1. 

It follows from l|26|) the phase and group velocities are determined from 

J£ = £ = cV-sinh^v 7 ! - ^. (28) 

Using the identity ismO — sinh(i6>) l|28|l becomes 

duj ui 



csin(^T^I), (29) 



where fi > 1. The latter can be approximated for z « 1 as czyfi — 1. In 
the general cosmos where fi > 1 gravity wave propagate with the velocity 
csin(^\/fi ~ 1) ~~ * c where v — > c. In the cosmos where fi < 1 we have evanescent 
decay. 

The dependence in l|16f) is not the same within a bound galaxy of stars 
and gas as it is for the large scale structure of the expanding universe, which 
considers only the center of mass motion of galaxies within it. This is because 
within a galaxy (or cluster) the full effect of the Hubble expansion is not felt with 
respect to the center of mass, of which the gravitational radiation must travel 
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with respect to, and here we show it results in the same solution as standard 
General Relativity. 

For spherically symmetric distribution of matter in a galaxy with the region 
of interest far from the central potential of a fixed mass, in the disk region, it 
may be derived (Eqs B.63a and B.67 of Carmeli 0) that 

where z = v/c and Q is the mass density, with the origin of coordinates coin- 
ciding with the origin of the spherically symmetric gravitational potential. 
Using (|30|l in (|15|l with (|13|l results in a wave equation 

( L^L — -— n ~ 1 f 2 92 — 

\c 2 dt 2 dr 2 r dr c 2 t 2 \ dr 2 dr 

where we can neglect the terms in the curly brackets because they are insignif- 
icant on the scale of a galaxy. Also after neglecting the third term for distant 
sources, we get the normal gravity wave equation of GR, that is, 

fid 2 d 2 \ - 

(?5P- 5?) < 32 » 

2.3 Density scales in the universe 

On the local scale the Hubble law v — r/r does not apply or is so insignificant 
as to be negligible. On that scale therefore effects of spacevelocity are negligible, 
or in other words, dv — ► 0. That is the realm where GGR reduces to the usual 
special and general relativity theory. Gravity waves propagate as is usually 
expected according to 

On the cosmological scale we expect the Hubble law to be very significant 
and hence that is the realm where spacevelocity is operative. On that scale we 
can now also consider 1)23(1 to represent a genuine modification to the usual 4D 
spacetime equation found in general relativity textbooks. The effect of spacevel- 
ocity is contained in the modified D'Alembertian operator. This means close to 
the source gravitational energy is emitted in the usual fashion as described by 
(|10|l . Over cosmological length scales, however, CGR predicts the gravitational 
waves from distant galaxies will be fully attenuated by the time the reach the 
Earth. Thus, projects such as the Large-scale Cryogenic Gravitational-wave 
Telescope (LCGT) in Japan could be very important to test CGR theory |12j . 
For example this project will detect gravity-waves from coalescing neutron-star 
binary systems 200 Mpc away at a SNR of 10. In contrast the TAMA field of 
view is 1 Mpc and LIGO is 20 Mpc. Thus, if the LGCT field of view includes a 
much larger volume of the Universe than LIGO or TAMA, and a smaller event 
rate than predicted by GR, this could be evidence for CGR. 

The best estimates of the local baryonic matter density puts it around Q m rs 
0.02 HU for the present epoch. As has been shown using the Carmeli 

theory it is not necessary to assume any dark matter in the cosmos, therefore 
the matter density between galaxies should be Q < 1 even out to a redshift 
z = 2 9 . Equation (|16|) is valid for large z therefore the analysis applies. 

Accordingly the averaged matter density of the universe has been determined 
at the current epoch Sl m = 0.021 ± 0.042 jS] and can be approximately related 
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by O = f2 m (l + z) 3 as a function of rcdshift for z < 1. Using this relation 
and the form of l|26l) the decay constant (k) is shown as ck/oj in Fig. ^ as a 
function of redshift, z. The value of ck/uj in Fig. ^is limited by the unknown 
form of fi(z) for z ^> 1. However, over the epochs shown, CGR predicts that 
gravitational waves do not propagate at scales beyond galaxies (and clusters). 

Gravity waves that are generated within a galaxy quickly decay in the void 
between. Gravitational radiation therefore leaks into the surrounding space 
according to l|23[) with attenuated amplitudes when Q drops below unity. Ac- 
cording to CGR, gravity waves will not propagate far in an expanding universe. 
Therefore we would expect to see no stochastic gravity wave background spec- 
trum. Instead we conjecture that the energy is deposited into space as heat. As 
a result they may contribute to the CMB blackbody temperature. 

The binary pulsar PSR B1913+16, discovered in 1974 by Russell Hulse and 
Joseph Taylor 11 , for which they won the 1993 Nobel prize, consists of two 
neutron stars closely orbiting their common center of mass. One of them is 
a pulsar with a rotational period of 59 ms and extremely stable compared to 
other pulsars. The two neutron stars slowly spiral toward their common center of 
mass radiating energy. The orbit period however is declining by about 7.5 x 10 -5 
seconds per year on an orbit period of 7.75 hours [T^]. This change, believed to 
result from the system emitting energy in the form of gravitational waves, has 
been a very precise and successful test of general relativity ^Ij • 

On the scale of the Galaxy it expected that there are still some small effects 
of the Hubble law EH ■ These effects modify the dynamics of the motion 
of tracer gases in the outlying regions of the Galaxy but in regard to gravity 
waves in galaxies they propagate unhindered, as in normal GR. Therefore it is 
expected that the energy dissipated by the Hulse- Taylor binary does travel as 
gravity waves within the Galaxy, but will be attenuated outside the region of 
the Galaxy where the mass density f2 drops below unity. 

3 Conclusion 

This paper derives the wave propagation equation for gravitational radiation in 
an expanding universe where an additional constraint has been placed on the 
nature of space itself. This is the introduction into the metric the fundamental 
assumption of the expansion of space according to the Hubble law. It is then 
found that no unattenuated gravity wave propagation may be possible in regimes 
where the Hubble expansion has effect. Then, depending on the density of 
matter, the propagation constant of any gravitational wave is either real or 
imaginary. If imaginary it represents an evanescent wave, which we conjecture 
means the energy is dissipated into the surrounding space as heat. When gravity 
waves are eventually detected, a test of this theory would be the detection of 
gravity waves from within the Galaxy but not from extra-galactic sources. 
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